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ABSTRACT 

Context. Ground-based high-dispersion (R~100,000) spectroscopy (HDS) is proving to be a powerful technique with which to char¬ 
acterize extrasolar planets. The planet signal is distilled from the bright starlight, combining spectral and time-differential filtering 
techniques. In parallel, high-contrast imaging (HCI) is developing rapidly, aimed at spatially separating the planet from the star. 
While HDS is limited by the overwhelming noise from the host star, HCI is limited by residual quasi-static speckles. Both techniques 
currently reach planet-star contrast limits down to ~ 10~ 5 , albeit for very different types of planetary systems. 

Aims. In this work, we discuss a way to combine HDS and HCI (HDS+HCI). For a planet located at a resolvable angular distance from 
its host star, the starlight can be reduced up to several orders of magnitude using adaptive optics and/or coronography. In addition, 
the remaining starlight can be filtered out using high-dispersion spectroscopy, utilizing the significantly different (or Doppler shifted) 
high-dispersion spectra of the planet and star. In this way, HDS+HCI can in principle reach contrast limits of ~ 1 CL 5 x 1 (L 5 , although in 
practice this will be limited by photon noise and/or sky-background. In contrast to current direct imaging techniques, such as Angular 
Differential Imaging and Spectral Differential Imaging, it will work well at small working angles and is much less sensitive to speckle 
noise. For the discovery of previously unknown planets HDS+HCI requires a high-contrast adaptive optics system combined with 
a high-dispersion R~ 100,000 integral field spectrograph (IFS). This combination currently does not exist, but is planned for the 
European Extremely Large Telescope. 

Methods. We present simulations of HDS+HCI observations with the E-ELT, both probing thermal emission from a planet at infrared 
wavelengths, and starlight reflected off a planet atmosphere at optical wavelengths. For the infrared simulations we use the baseline 
parameters of the E-ELT and METIS instrument, with the latter combining extreme adaptive optics with an R= 100,000 IFS. We 
include realistic models of the adaptive optics performance and atmospheric transmission and emission. For the optical simulation we 
also assume R= 100,000 IFS with adaptive optics capabilities at the E-ELT. 

Results. One night of HDS+HCI observations with the E-ELT at 4.8 pm (A A = 0.07/jm) can detect a planet orbiting orCen A with a 
radius of R=1.5 Re^th and a twin-Earth thermal spectrum of T cq =300 K at a signal-to-noise (S/N) of 5. In the optical, with a Strehl 
ratio performance of 0.3, reflected light from an Earth-size planet in the habitable zone of Proxima Centauri can be detected at a S/N 
of 10 in the same time frame. Recently, first HDS+HCI observations have shown the potential of this technique by determining the 
spin-rotation of the young massive exoplanet /? Pictoris b. 

Conclusions. The exploration of the planetary systems of our neighbor stars is of great scientific and philosophical value. The 
HDS+HCI technique has the potential to detect and characterize temperate rocky planets in their habitable zones. Exoplanet scientists 
should not shy away from claiming a significant fraction of the future ELTs to make such observations possible. 

Key words. Infrared: planetary systems - Methods: data analysis - Techniques: imaging spectroscopy - Techniques: high angular 
resolution 


1. Introduction 


Since the late 1990s it has been recognized that high-dispersion 
spectroscopy (HDS) could be a powerful way to characterize 
extrasolar planet atmospheres. First attempts to use HDS were 
focussed on detecting optical starlight reflected off hot Jupiter 
atmospheres. However, their optical albedos turned out to be 
low (e.g. Charbonneau et al.|1999[ Collier-Cameron et al.|2004| 


Leigh et al.|2003af Leigh et al.|2003b)>, making these measure- 


with HDS also resulted in upper limits (Wiedemann et al. 2001 
|Deming et al.|2005] ). 

Recently, the availability of more stable infrared spectro¬ 
graphs, in particular the CRyogenic high-resolution InfraRed 
Echelle Spectrograph (CRIRES; |Kaeufl et al.||2004| on the 
Very Large Telescope (VLT, European Southern Observatory), 
in combination with new observational strategies, has led to the 
first detections. |Snellen et al.| ( [2010| ) measured carbon monoxide 
at 2.3 pm in the transmission spectrum of HD 209458b, reveal- 


ments very challenging. Early attempts to detect molecular ab¬ 
sorption features in the thermal emission spectra of hot Jupiters 


ing the orbital motion of the planet. Brogi et al. (2012 1 found 
the same molecular species in the thermal spectrum of the non- 
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transiting planet r Bootis b, determining its mass and orbital in¬ 
clination (see also Rodler et al.|201 2), and it was also detected in 


the thermal spectra of HD 189733 b (|de Kok etal.|2013[ Rodler 


et al. 2013|> an d tentatively 51 Pegb (Brogi e t al,|2013[ >. Recently, 
Birkby et al. (|2013j> measured water absorption in the thermal 


spectrum of hot Jupiter HD 189733 b at 3.2 pm using HDS with 
CRIRES, and jLockwood et ah| ( [2014j ) in that of r Bootis b using 
the Keck telescope. 

In parallel to HDS, high-contrast imaging (HCI) is rapidly 
developing, driven by improvements in adaptive optics technol¬ 
ogy and coronographic concepts. This technique is now well es¬ 
tablished for young, massive self-luminous planets at significant 


orbital distances from their host star (e.g. 

Chauvin et al. 20051 

Lagrange et al. 2009. 

Marois et al. 2008 

Kalas et al. 2008). 

Several new HCI facilities are coming on 

ine, such as GPI on 

the Gemini Telescope (Macintosh et al. 2014), SPHERE on the 


Very Large Telescope (jBeuzit et al. 2Q10j ), and SCExAO on Sub¬ 
aru (Jovanovic et a l,|2013| >,~ which will push these observations 
to closer and cooler planets. 

In this paper we discuss a way to combine high dispersion 
spectroscopy with high-contrast imaging (HCI). For a planet 
located at a resolvable angular distance from its host star, the 
starlight can be reduced up to several orders of magnitude us¬ 
ing adaptive optics and/or coronography. In addition, the re¬ 
maining starlight can be filtered out using high-dispersion spec¬ 
troscopy, utilizing the significantly different (or Doppler shifted) 
high-dispersion spectra of the planet and star. In Sect. 2 we give 
technical descriptions for HDS, HCI, and HDS+HCI. In Sect. 3 
we describe the instrument requirements for HDS+HCI, which 
basically come down to a high-contrast adaptive optics system 
combined with a high-dispersion R~ 100,000 integral field spec¬ 
trograph (IFS). This combination currently does not exist, but is 
planned for the European Extremely Large Telescope. In Sect. 
4 we describe simulations of infrared HDS+HCI observations 
with the METIS instmment on the European Extremely Large 
Telescope (E-ELT), and in Sect. 5 that of optical HDS+HCI ob¬ 
servations with a hypothetical instrument on the same telescope. 
The results are discussed in Sect. 6. 

2. Technical descriptions 

2.1. High dispersion spectroscopy (HDS) 

High dispersion spectroscopy makes use of the fact that at high 
spectral dispersion molecular bands resolve into tens to hundreds 
of individual lines, which shift in wavelength during observa¬ 
tions thanks to a change in the radial component of the orbital 
motion of the planet. For hot Jupiters, which have orbital veloc¬ 
ities of >100 km sec -1 , this change in wavelength can be up to 
a several times the spectral resolution per hour. This is used to 
separate the planet spectrum from the quasi-stationary telluric 
and stellar components which strongly dominate the observed 
spectra. The measurement of a molecular species is subsequently 
enhanced by cross-correlating the observed spectra with a model 
template spectrum, such that the signals from the individual lines 
are optimally combined (e.g. Brogi et al.|2013| >. 

Neglecting the interference by telluric absorption, the 
achieved signal-to-noise ratio (S/N) of HDS is to first order given 
by 

S/N = , Sp ‘ anet — VAW (1) 


4 


s 

star + °T g + °T 


4 - rr 

RN Dark 


where S p | ;lncl is the planet signal, S star is the signal from the star 
(both in units of photons per resolution element), and <r bg , <xrn. 
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and cr Dalk are the photon shot noise from the sky and telescope 
background, the readout noise, and the noise from the Dark cur¬ 
rent respectively. Mines is a multiplication factor that takes into 
account the number and strength of the individual planet lines 
targeted. For most HDS observations, <x bg , <x RN and cr Dark can 
be neglected because the stars are bright, and the backgrounds 
low. Combining the thousands of lines in the optical spectrum 
of a sunlike star typically results in a S/N enhancement of a fac¬ 
tor of ~30 (Mii n es ~ 1000, |Collier-Cameron et al.|2004| l. In this 
way, |Leigh et al.[ ( f2003a| ) reached a lcr planet/star contrast limit 
of 1.3x1 O' 0 for the HD75289b system (V=6.3) by combining 
684 spectra taken over 4 nights with UVES on the VLT. For r 
Bootis b (K=3.36) at 2.3 pm. [Brogi et al.| ( j2012) reached a lcr 
CO line-contrast of ~2xl0 s with respect to the stellar spectrum 
in 15 hours of CRIRES VLT observations, utilizing a factor of 
~5 enhancement by combining the signal from two-dozen CO 
lines. [Birkby et ak| ( |2013| l reached a lcr H 2 O line contrast of 
<3x10 -4 in 5 hours of observations at 3.2pm for HD189733 b 
(K=5.54) using the same instrument. Because of the strong tel¬ 
luric absorption in this wavelength regime, the S/N enhancement 
from combining numerous but weak water lines is only a factor 
of ~3. _ 

[de Kok et al.[ ( |2014) > identified new opportunities of using 
CRIRES to probe hot Jupiters, showing that several other wave¬ 
length regions promise to give strong molecular signals, such 
as around 3.5 pm for day side spectroscopy and at 1.6 pm for 
transmission spectroscopy. In addition, it is expected that in 
certain cases the planet night-side may even produce stronger 
signals than its dayside. Although the nightside temperature is 
lower and therefore the planet signal fainter, the absorption lines 
themselves may be stronger as a result of a steeper atmospheric 
temperature-pressure profile. 

Snellen et al. ( |2013| l explored future possibilities of high- 
dispersion transmission spectroscopy to detect molecular oxy¬ 


gen in nearby Earth-like exoplanets (see also Rodler & Lopez- 
|Morales|2013[ ). They showed that the transmission signal of the 
oxygen A-band from an Earth-twin orbiting a small red dwarf 
star is only a factor of 3 smaller than that of carbon monoxide 
detected by |Brogi et ~aL] ( |2Q12) > in the hot Jupiter tau Bootis b, 
albeit such a star will be orders of magnitude fainter. If Earth¬ 
like planets are very common, the planned extremely large tele¬ 
scopes could detect an oxygen signal within a few dozen transits, 
but this could take decade(s), since only a few transits would be 
favorably observable from one location on Earth per year. De¬ 
tection of an oxygen signal from an Earth-like planet in the hab¬ 
itable zone of a solar type star could take centuries, since such 
transits are so rare. We note that this same idea, HDS of the oxy¬ 
gen A band, was already suggested more than a decade ago by 
|Webb & Wormleaton| ( |200T) , but they argued it could take 10 
hours or less on an 8m telescope. 

2.2. High-contrast imaging (HCI) 

Planet detections with HCI lead to planet characterization when 
observed at multiple wavelenths. Broadband photometry have 
lead to estimates of planet effective temperatures and radii, and 

Currie 


possible prevalance of clouds (e.g. Currie et al. 2011 


|et al.|[20l3[ |Bonnefoy et al.||2013] >. However, for ground based 
observations, the turbulent atmosphere degrades the diffraction 
limited images of stars into a seeing limited disk. All large tele¬ 
scopes use an adaptive optics system to measure the atmospheric 
distortion with a wavefront sensing (WFS) camera and then ap¬ 
ply optical correction with a deformable mirror to drive the 
sensed wavefront aberrations to a null. The resultant diffracted 
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halo of light from the primary star dominates the noise contribu¬ 
tion at the location of the exoplanet (see middle panel of Figure 
1). In addition, the optical WFS camera path aberrations and sci¬ 
ence instrument optical path aberrations vary slowly with time, 
changing the relative intensity of the diffraction halo at a level 
of 10 4 to 10 5 of the stellar peak flux. Differencing two im¬ 
ages taken several minutes apart shows the presence of quasi¬ 
static speckles similar in angular size to potential faint compan¬ 
ions. These are attributed to the time-varying nature of the non¬ 
common path aberrations between the adaptive optics system 
WFS camera and science camera. Several high contrast observ¬ 
ing techniques and post-processing algorithms have been devel¬ 
oped to form a two-step process to estimate the science PSF as 
a function of time during the science exposures that gain extra 
orders of magnitude in planet/star contrast. These include Angu¬ 
lar Differential Imaging (ADI; Marois et al.|2006| ), Spectral Dif¬ 
ferential Imaging (SDI; Rac ine et al.|1999), Locally Optimized 
Combinations of Images (LOCI; Lafreniere et al. 2007), and 


Principal Component Analysis (PCA; Amara & Quanz 


2012] 


|Soummer et al.|2012[ |Meshkat et al.|2014) . Some of these tech¬ 

niques can also be used in combination to get the best contrast. 

Coronagraphs are angular filters that reject on-axis light 
whilst transmitting light from nearby off-axis faint companions, 
reducing the diffraction halo from the primary star through op¬ 
tics introduced in the focal plane and pupil planes of the science 
camera (see Guyon et al. 2006a| for a review). Coronagraphs re¬ 
duce the flux in diffraction structures and so the effects of non¬ 
common path errors are correspondingly reduced. Several coro- 
nagraphic designs for space-based high contrast designs (le-9 to 
le-10) have been developed (Phase-induced amplitude apodiza- 
tion or PIAA - Guyon et al.|2006b) Vortex Coronaraph, Mawet 


et al j |2009| Band Limited Coronagraph, Vanderbei et al. |2003[ 
Cash 2006) but all are sensitive to pointing errors and vibrations 
that the AO systems cannot remove, and are chromatically lim¬ 
ited. For ground based telescopes, the raw constrast is limited to 
ICC 5 at a few ,1/D because of the time delay between the WFS 
and the AO system response, (Guyon et al.|2012). In these cases, 


pupil plane coronagraphs that use phase apdization (APP; 


wort hy et al.|2007) |Codona & Angel|2004| broadband APP Ot- 


Ken- 


ten et al. 2014| l or amplitude apodization (Carlotti 201 3 ) can pro- 
vide robust suppression at small angular separations ( |Kenwo rthy 
|etal.[2 013; Meshka fet al.|2014] > in science instruments without a 
dedicated coronagraphic optical design. An advanced APP with 
a limited dark area in the shape of a slot can be designed with a 
contrast of 10 10 or better with a transmission of 95% (Keller et 
al. in preparation). When manufactured as a vector APP (Snik et 
al. 2012), such coronagraph will work over a wavelength range 
of an octave while being largely immune to any remaining tip-tilt 
errors in the incoming wavefront. 


2.3. Combining HDS and HCI 


|Sparks & Ford| (2002 ) were the first to advocate a method that 
combines coronagraphic imaging with what they called “spec¬ 
tral deconvolution”, using an integral field spectrograph (IFS; 
see also [Riaud & Schneider] ( j2007[ )). The principles laid out in 
this pioneering work are very similar to what we present here. 
The technique introduced by Sparks & Ford (2002) has recently 
been implemented by |Konopacky et al. ( 2013) 1, who have de¬ 
tected carbon monoxide and water in the spectrum of the directly 
imaged exoplanet HR8799c. An important aspect of the method 
described below is the high spectral dispersion. This makes it 
possible to distinguish between direct starlight and Doppler- 
shifted starlight reflected off a planet atmosphere. It can also be 


used to distinguish between the telluric aborption spectrum and 
Doppler-shifted absorption lines in the thermal spectra of plan¬ 
ets very similar to our Earth. Kawahara et al.]( |2014| ) propose a 
related technique, combining coronagraphy with HDS to probe 
hot close-in planets such as hot super-Earths. 


For a planet located at some angular separation from its 
host star, the starlight at the position of the planet can be re¬ 
duced up to several orders of magnitude using adaptive optics 
and coronography. Subsequently, the remaining starlight can be 
filtered out using high-dispersion spectroscopy, utilizing the sig¬ 
nificantly different (or Doppler shifted) high-dispersion spectra 
of the planet and star. While in classical HDS it is the change 
in the Doppler-shift of the planet signal that is used to separate 
it from the much stronger (quasi-)stationary stellar and telluric 
contributions, in the case of HDS+HCI it is the fact that the star- 
spectrum and telluric absorption are dominant but identical over 
the whole field, while the planet signal itself is uniquely different 
and strongly localized at one position in the image. The speckles 
have the high-dispersion spectrum of the star and are therefore 
effectively filtered out, influencing only the photon-noise statis¬ 
tics. In contrast to image processing techniques such as ADI and 
SDI, HDS+HCI will work well at small working angles (a few 
d/D) and is much less sensitive to quasi-static speckle noise that 
tends to plague HCI. 


While both HDS and HCI can each reach (raw) planet/star 
contrasts down to ~ 10~ 5 , a combination of HDS+HCI can in 
principle reach contrasts of ~ ICC 5 x ICC 5 , although in practice 
this will be limited by photon statistics and/or sky-background. 
We also note that at a dispersion of R= 100,000 the scattered 
starlight may be coherent, leading to high-frequency fringes, 
i.e. the speckle pattern as a function of wavelength may not be 
smooth. This will need to be tested to assure that the above as¬ 
sumptions are correct. 


The principle of HDS+HCI is illustrated in Fig. 1, showing 
in the left panel a toy-model of the stellar point-spread function 
(PSF) for seeing-limited observations (seeing = 0.6 arcseconds). 
A hot Jupiter currently can be detected at a contrast of < 10 4 in 
seeing-limited observations, a level of contrast that has already 
been achieved both in the optical (e.g. Leigh et al.|2003a I and in¬ 
frared (e.g. Brogi et al.|2012) . The middle panel of Fig. 1 shows 
the PSF of AO-assisted HCI observations with an 8m telescope 
at 0.5/rm, with a Strehl ratio of 0.3 under 0.6 arcsecond seeing 
conditions. In this example, a raw contrast of 1 x 10 5 can be 
achieved 0.6 arcseconds away from the host star (without utiliz¬ 
ing data analysis techniques such as SDI or ADI). This means 
that by combining HDS+HCI, a contrast of 10 3 x 10 4 = 10 7 
can be achieved at the planet position, as shown in the right panel 
of Fig. 1. In addition, each detection of a planet automatically 
results in a measurement of the radial component of its orbital 
velocity. By monitoring a system over time, both the planet posi¬ 
tion and the change in its radial velocity can be used to determine 
its orbital parameters, including the inclination of the system. 


The contrasts quoted above provide an indication of the 
currently proven limits of HDS+HCI before systematic effects 
could start to play a role because both such raw imaging con¬ 
trasts and HDS signals at these low levels have already been 
reached. With the starlight reduced by a factor K at the planet po¬ 
sition (compared to when no attempt is made to angularly sepa¬ 
rate the planet from the star), the S/N reached within a certain ex¬ 
posure time will be reduced by a factor yK, and for HDS+HCI 
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High Dispersion Spectroscopy (HDS) High—controst Imaging (HCI) HDS + HCI 





Fig. 1 . Toy model of the HDS+HCI method. The left panel shows the stellar point-spread function (PSF) for conventional, seeing-limited (see- 
ing=0.6 arcseconds) HDS observations, with indicated a contrast of lxlO -4 . This level of contrast has been readily achieved both in the optical 
(e.g. |Leigh et al.|2003aj l and infrared (e.g. Brogi et al. 2012[ >, meaning that such planet signal can be detected at a 1:10,000 level in the spectrum 
of the star. The middle panel shows a model PSF for HCI observations for an adaptive-optics assisted 8m telescope with a Strehl ratio of 0.3 at 0.5 
pm, under the same seeing conditions. The PSF is modeled as the theoretical Airy profile of the telescope combined with a Moffat function as the 
non-AO-corrected seeing-limited contribution. A hypothetical planet is inserted at an angular distance of 0.6 arcseconds from the star at a contrast 
of 1:1,000 with respect to the stellar brightness at the position of the planet. The right panel illustrates that in this example HDS+HCI can achieve 
a contrast of 10 -3 x 10 4 = 10 -7 at the planet position. 


Eq. 1 becomes 

S/N = Splanet VAW^, (2) 

JS stal /K + crl s + o-l N + 0"Dark 

meaning that for the same telescope and exposure time, the 
planet/star contrast improves by a factor \[K with respect to 
classical HDS. In addition, the higher K, the higher the relative 
importance of the sky background (<Xb g ), which in particular be¬ 
yond 3 pm this becomes the dominant noise factor in Eq. 2. It 
could mean that in the background limited regime existing HCI 
techniques will be more sensitive than HDS+HCI, depening on 
the strength and number of planet absorption lines targeted in 
the observed wavelength regime, and on the relative instrument 
through-put and bandwidth. In any case, HDS+HCI will have as 
added value that it also measures the radial component of the 
orbital velocity of the planet, helping to determine its orbital pa¬ 
rameters. 

Imaging contrast will be pushed significantly further with 
the next generation of Extremely Large Telescopes (ELTs). At 
the diffraction limit, angular resolution scales with A/D, where 
D is the telescope diameter, meaning that the angular distance 
in units of diffraction elements increases linearly with D. The 
relevant parameter for these exoplanet observations is the con¬ 
trast of a point source compared to the local intensities of the 
star and background. The FWHM of the central part of the Airy 
disk scales with D -2 , while the relative surface brightness of the 
uncorrected (seeing-limited) part of the PSF is invariant. This 
means that for the same Strehl ratio, the stellar halo (and the 
sky background) are reduced by D 1 2 , while the collecting power 
scales as D 2 . In other words, in the context of a Nyquist-sampled 
PSF at the diffraction limit of the telescope, the increase in diam¬ 
eter implies that two effects are combined; 1) a modification of 
the PSF shape for a given Strehl leading to a flux concentration 
in a PSF core scaled with D 2 , and 2) a collecting power scaled 
also as D 2 . Furthermore, the influence from the star is likely to 
decrease for larger D because the planet-star separation is larger 
in terms of A/D. Altogether, a certain planet/star contrast can be 
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reached in a time that scales as 1/D 4 or more with the telescope 
size. 


3. Instrument requirements for HDS+HCI 


Assuming that the positions of the planets are unknown, 
HDS+HCI requires an adaptive-optics assisted telescope with 
an Integral Field Spectrograph (IFS). The IFS could either be 
fiber-fed or use image slicing techniques. In the search for re¬ 
flected starlight, the optimal spectral dispersion is set by the in¬ 
trinsic width of the stellar lines, which for slow rotators is dom¬ 
inated by the turbulent velocity fields in the stellar photosphere 
of ~ 4 km sec -1 . This provides an optimal spectral dispersion of 
R=75,000 in the optical. The intrinsic widths of absorption lines 
in a planet thermal spectrum can be significantly more narrow, 
but integrated over the planet disk are likely to be broadened by 
the planet’s spin velocity. All rocky planets in the Solar System 
have spin velocities of <1 km s -1 , while those of the gas giants 
extend up to 12.5 km s -1 for Jupiter (and 25 km s -1 for f3 Pic- 
toris b; Snellen et al. 2Q14| >. Experience with CRIRES at the VLT 
shows that a high dispersion is also important for the effective 
removal of telluric contaminations. It indicates that R= 100,000- 
150,000 is a good compromise. 

The larger the instantaneous spectral coverage, the higher the 
enhancement in S/N will be from combining the signals from in¬ 
dividual absorption lines, as defined by /Vi mcs in Eq. 1 & 2. For 
probing optical reflected starlight, the 0.4-0.7 pm wavelength 
regime, as commonly used by modern optical echelle spectro- 


1 Already such an instrument is planned, called METIS - the Mid- 
infrared ELT Imager and Spectrograph for the E-ELT ( jBrandl et al.| 
|201 Of Its baseline design includes a high-dispersion R= 100,000 inte¬ 
gral field unit for L (2.9-4.2 pm) and M (4.5-5.0 pm) band spectroscopy 
with an instantaneous wavelength coverage of 0.07pm, and a field of 
view of 0.4xl.5 arcseconds cut in 24 slices. To our knowledge, this is 
the first (planned) instrument that includes a high-dispersion IFS, on a 
telescope that will utilize adaptive optics techniques. It will be ideally 

suited for HDS+HCI observations. It may also be equipped with a long- 
slit unit delivering an instantaneous wavelength coverage of ~0.5 pm - 
a powerful addition for planets with a known position. 
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Table 1. Main assumptions for the infrared HDS+HCI simulations as 
presented in Figs. 2 and 3. 


Planet Spectrum 


Telescope + Instrument 
Telescope collecting area 
Telescope temperature 
Telescope emissivity 
Telescope+instrument throughput 
AO Strehl (4.85 fim) 

Spectral resolution 
Exposure time 
Spectral range 

Target: a Cen A 
Apparent K magnitude 
T e ff (star) 

Stellar radius 
Distance 
Planet radius 
Planet radial velocity 
T e ff (planet) 

Bond albedo 
Planet spectrum 


976.3 m 2 
280 K 
0.15 
15% 
0.9 

R= 100,000 
30 hr 
4.82-4.89 pm 


-1.47 
5800 K 
1.22 R sun 
1.34 pc 
1.5 R Earth 
30 km sec -1 
300 K 
0.3 

Earth-like 



Stellar Spectrum + Telluric absorption 



4.83 4.84 4.85 4.86 4.87 4.88 

Wavelength (jj. m) 


graphs, will cover the densest regions in solar-type spectra. Do 
note however that for an IFS with R=75,000 this will require 
a large number of detectors. For cooler host stars, the red part 
of the optical spectrum is more interesting, up to the 1pm CCD 
sensitivity cut-off, where most of the energy of M-dwarfs is emit¬ 
ted, and the molecular bands from titanium oxide and vanadium 
oxide are prominent. In addition, if Earth-like planets could be 
probed, the oxygen A-band at 0.76 pm would be extremely in¬ 
teresting. Similar arguments exist for targeting the intrinsic ther¬ 
mal spectra of exoplanets. Although the entire wavelength range 
up to 5 pm is interesting for ground-based HDS (beyond that 
the sky background is too high), practical limits in wavelength 
coverage are governed by detector technology and windows of 
transparency in the Earth atmosphere. The strongest signatures 
can be expected from water, methane, carbon dioxide, and also 
from carbon monoxide for the hottest planets. 

The field of view of the IFS should be at least on the order of 
~ 1 — 2 arcseconds to include the habitable zones of the nearest 
solar-type stars. Particularly interesting for optical HDS+HCI 
are the small inner work angles, meaning that for the ELTs the 
habitable zones of nearby M-dwarfs at 0.05-0.2 AU could be 
within reach. This science case allows a significantly smaller 
field of view. Another important consideration for the instru¬ 
ment design is that typically, image sheer IFUs are not isotropic, 
meaning that they provide better resolution along the slices than 
across the slice direction. Hence, one should observe half of the 
time with the IFU rotated 90 degrees, and optimally combine the 
two orientations in software. We note that both the contrast from 
reflected starlight and thermal emission are expected to decrease 
significantly as a function of orbital distance. Although young, 
self-luminous planets can be seen much farther out, after discov¬ 
ery their positions are known and can therefore also be probed 
with a long-slit HDS to determine their radial velocity and atmo¬ 
spheric constituencies. 


4. Simulations of infrared HDS+HCI observations 

We simulated infrared HDS+HCI observations with the E-ELT 
(see Table[I]) of hypothetical planets orbiting a Cen A (K=-1.47, 
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Fig. 2. Top panel: Simulated E-ELT planet spectrum as photons per 
Av = 3 km s -1 wavelength interval. The dotted line indicates the Earth 
telluric spectrum, which in this case is very similar to the simulated 
thermal spectrum of the exoplanet (solid line), with the latter shifted by 
30 km sec -1 in velocity because of its orbital motion around trCen A. 
Depending on the time of year, the telluric absorption may show an ad¬ 
ditional ±30 km sec -1 velocity with respect to the exoplanet because of 
the Earth’s orbital velocity around the Sun. Middle panel: The stellar 
spectrum in photons per wavelength step as seen through the Earth at¬ 
mosphere, at the sky position of the planet. All the absorption lines are 
telluric in nature. Bottom panel: The background spectrum in photons 
per wavelength step per pixel, modeled as contributions from the night 
sky and emissivity of the E-ELT telescope. The former consists for a 
large part of bright emission lines at the wavelengths of the telluric ab¬ 
sorption lines. 


T. st ar=5800 K, R star =1.22 R sun , dist= 1.34 pc). The planet radii 
were varied between 1 and 1.5 R^arth. with equilibrium tem¬ 
peratures between 255 K (like our Earth) and 355 K. The as¬ 
sumed Bond albedo of Ab=0.3 implies orbital radii in the range 
0.63-1.22 AU, corresponding to angular distances in the range 
0.48-0.93 arcsec. The planet’s orbital velocity was fixed to 30 
km sec -1 , similar to that of the Earth. The planet spectrum (Fig. 
3) was modeled as the thermal spectrum of the Earth scaled 
by the size of the planet disk, and includes molecular absorp¬ 
tion from carbon dioxide, methane, and water vapor. All atmo¬ 
spheric parameters, T/p profile and molecular abundances were 
fixed to the Earth value. No attempt was made to introduce atmo¬ 
spheric chemistry when changing the atmospheric temperature. 
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CC map - a Cen A with METIS 



0.4 



0.0 0.5 1.0 

Distance (arcsec) 


Fig. 3. Top panel: Slice at the V=30 km sec -1 of the cross-correlation 
data cube showing the simulated planet (R p = 1.5 R Eart h and T eq =300 
K) around a Cen A, for a 30 hours exposure with METIS at the E- 
ELT. The planet is visible 0.7 arcseconds away from the star at an S/N 
of 8. Since the cross-correlation technique automatically renormalizes 
the noise, the star is not visible in this data cube. The bottom panel 
shows the same but now multiplied by the local contrast to emulate a 
conventional HCI image. 


The Earth model spectrum was simply scaled according to the 
change in the blackbody emission at the observed wavelength. 

We used the baseline instrument parameters of METIS at the 
E-ELT ( |Brandl et al.|20l2| . The E-ELT was assumed to have a 
collecting area of 976.3 nr, at a temperature of 280 Kelvin with 
an emissivity of 0.15. The total telescope-tinstrument through¬ 
put (excluding atmosphere) was assumed to be 15%. We set the 
wavelength range to 4.82-4.89 pm, optimizing the combination 
of strong planet absorption lines and weak telluric absorption. 
The Earth atmospheric emission and transmission spectra were 
taken from a recent sky model for ESO exposure time calcu¬ 
lator^] The AO-assisted PSF for METIS has been modeled by 
Stuiket al. (2010). At this wavelength the extreme AO system on 
the E-ELT is expected to deliver a Strehl ratio of 0.9 for bright 
stars. A simulated point-spread function (PSF) for this level of 
correction, kindly provided by the METIS instrument team, was 
used for our simulations. The total exposure time varied from a 
few hours to several nights to determine the S/N limits of these 
observations. 

We first calculated the observed flux from the solar-type 
star in photons per wavelength interval, and multiplied it with 
an Earth telluric transmission spectrum. This was subsequently 
convolved with the model PSF to provide a data cube of stel- 


2 See http://www.eso.org/observing/etc/doc/skycalc/helpskycalc.html 
#version 
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Contrast curve 



Distance (arcsec) 

Fig. 4. The contrast curve for the simulated infrared observations, as 
presented in Figures 3 and 4, and described in section 4, for an obser¬ 
vation of a Cen A with METIS for 30 hours centered at a wavelength 
of A = 4.85yum. The asterisk signifies the simulated planet detected at 
an S/N of 8. The dotted lines show scaling laws for planets with the 
same spectrum, but different effective temperatures (T eq = 250, 300, 
and 350 K) and radii (R = 1.0 and 1.5 REarth)- No attempt is made to 
include changes in chemical composition or temperature structure for 
such planets. 


lar flux per spaxel (wavelength interval per IFS pixel). The same 
was done for the planet, by appropriately scaling down the planet 
spectrum and Doppler shifting it to the correct wavelength scale, 
and offsetting the planet position by the angular separation cor¬ 
responding to its equilibrium temperature. The sky background 
spectrum was scaled to the collecting area of the E-ELT, the 
pixel size, and wavelength intervals. This was added to the pho¬ 
ton flux expected from the emissivity of the telescope itself. In 
this way three data cubes were constructed, one with the pho¬ 
ton flux from the star, one with the flux from the planet, and one 
with the flux from the background. Both the background (uni¬ 
form over the detector) and the contribution from the star (same 
spectrum over the detector but varying in amplitude) were subse¬ 
quently assumed to be removed during the data analysis in such 
way that only Poisson noise is left from their contributions. We 
note that this is an idealistic case. Although it has been shown 
that the HDS method is capable of removing stellar contribu¬ 
tions down to < ICE 4 level, this has not been shown to work 
down to this level in combination with high contrast imaging, 
and much remains to be understood for the subtraction of stel¬ 
lar speckles in the high-spectral resolution regime. We used the 
Gaussian approximation for their noise contributions by com¬ 
puting the square root of these data cubes, which were subse¬ 
quently multiplied by a random Gaussian noise distribution and 
added to the planet data cube. The spectrum at each pixel po¬ 
sition is subsequently cross-correlated with the planet template 
spectrum to obtain a three-dimensional data cube, containing x,y 
sky position, and the cross-correlation signal for velocities be¬ 
tween -150< V < +150 km sec' 1 . 

Figure 3 shows a slice of the cross-correlation data cube at 
the planet velocity (top panel) and a renormalized map (emu¬ 
lating a conventional HCI image - lower panel) for simulations 
of a planet with R p = 1.5 REarth and T eq =300 K for 30 hours of 
observing time. The planet is clearly visible -0.7 arcsec away 
from the star in x-direction, and is detected at an S/N of -8. 
In one night (10 hours) it would be detected at an S/N of -5. 
Figure 5 shows the contrast achieved as a function of angular 
distance, with the planet (R p = 1.5 REarth and T eq =300 K) of Fig. 
4 indicated by the asterisk. The dotted lines indicate scaling re- 
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Fig. 5. Graphical representations of the cross-correlation data cube for 
METIS observations. All parameters were kept the same as for the 
simulation presented in section 4. However, for illustrative purposes, 
the exposure time was increased such to show the planet signal at a 
S/N of ~20, and the planet velocity was increased by a factor of 2 to 
enhance the offset in velocity space. The cross-correlation data cube 
is renormalized as in the bottom panel of Figure 3. A short movie 
from this data cube, rotating in sky and velocity angle, is available at 
http ://www. strw.leidenuniv.nl/~ snellen/ 


lations for these simulations, for planets with radii of 1.0 and 1.5 
REarth, and equilibrium temperatures of T eq =250, 300, and 350 
K. Again, these assume the same planet spectra, only scaled in 
flux by the brightness temperature and planet size, and do not 
take into account expected changes of chemical composition or 
atmospheric structure as a function of temperature, and should 
therefore only be considered as an approximation. An Earth-like 
planet with T eq =250 and R p = 1.5 R e has a S/N of ~ 1 in these 30 
hour observations. 


CC map — Proxima 
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Fig. 6. HDS+HCI cross-correlation map of 10 hours of optical obser¬ 
vations with the E-ELT using a R= 100,000 IFS and an adaptive optics 
system producing a Strehl ratio of 0.3. The hypothetical planet with a 
radius of R=1.5 R Eart h, albedo of 0.3, and T eq =280 K such that it is at 
an orbital radius of 0.032 AU, 25 milliarcseconds from the star. The 
starlight reflected off the planet is detected at an S/N of ~10. 


Figure 5 shows a graphical representation of the cross corre¬ 
lation data cube for the simulated METIS observations. All pa¬ 
rameters are the same as above, except that for illustration pur¬ 
poses the exposure time was increased to more clearly show the 
planet (now at an S/N of ~20), and the planet velocity was dou¬ 
bled to enhance the offset in velocity space. 

If we assume that the targeted planet has already been dis¬ 
covered, e.g. using classical HCI techniques, a long-slit mode 
can be used to probe the spectrum of the object. This has the 
potential advantage of a significantly large instantaneous wave¬ 
length coverage, e.g. of 0.5 p m, instead of 0.07 pm for METIS. 
We found that this increases the S/N for a given exposure time by 
about a factor of ~2. We note however that it may be challeng¬ 
ing to center a planet, even when previously imaged, on a 10-20 
milliarcsecond wide slit of a single-slit spectrograph (a typical 
planet around a nearby star could move by this angular distance 
in one week). The utilization of slit-view cameras or automated 
centering only works for bright sources, and one would not want 
to integrate for many hours missing the planet. This issue does 
not exist for an IFS. 


5. Simulations of optical HDS+HCI observations 

We also simulated optical HDS+HCI observations with the E- 
ELT, for starlight reflected off planet atmospheres. Although 
currently the optical high-contrast imaging instrument for the 
E-ELT, EPICS ( |Kasper et al. |2010[ >, is proposed with a low- 
dispersion IFS, we would like to make a case for a high- 
dispersion IFS. At these wavelengths the angular resolving 
power of the E-ELT will be well suited to target planets in the 
habitable zones of nearby M-dwarfs. It is clear that an extreme 
adaptive optics system for the optical wavelength regime will be 
technically challenging ( Kasper et al.|2010| , and we therefore 
assume a lower Strehl ratio of 0.3 under 0.6 arcsec seeing con¬ 
ditions. For the integral field spectrograph we assume a spec- 

Article number, page 7 of[9] 










A&A proofs: manuscript no. 3M_snellen_accepted 


Table 2. Main assumptions for the optical HDS+HCI simulations as 
presented in Figs. 6. 


Telescope + Instrument 

Telescope collecting area 976.3 m 2 

Telescope+instrument throughput 15% 

AO Strehl (0.75 pm) 0.3 

Spectral resolution R= 100,000 

Exposure time 10 hr 

Spectral range 0.6-0.9 pm 

IFU pixels 30x30 2 mas 


Target: Proxima Cen 

Apparent V magnitude -11.05 

T eff (star) 3040 K 

Stellar radius 0.141 R sun 

Distance 1.30 pc 

Planet radius 1.5 REarth 

Planet radial velocity 30 km sec 1 

T e ff (planet) 280 K 

Grey geometric albedo 0.3 

Orbital radius 0.032 AU 

Angular distance from star 25 mas 


tral dispersion of R=100,000 with 30x30 elements of 2.0 mil- 
liarcsecond pixels. Hence the field of view is 60x60 mas, well 
matched to resolving the habitable zones of nearby M-dwarfs. 
The wavelength range is set to 0.6-0.9 pm. Extending this to 
shorter wavelengths does not significantly increase the sensitiv¬ 
ity because of the red colors of the host stars. Also, we note that 
this set-up will require several tens of CCDs. 

In our simulations we target our nearest neighbor Proxima 
Centauri (V=11.05 mag, T star =3040 K, R star =0.141 R sun , Dist= 
1.30 pc). We assume a hypothetical planet with a radius of R=1.5 
REarth, a gray albedo of 0.3, and T eq =280 K such that it is at an 
orbital radius of 0.032 AU, 25 milliarcseconds from the star. We 
assume that half of the planet is seen to be illuminated by its host 
star, as expected at the greatest elongation. 

We subsequently follow a similar strategy as for the METIS 
infrared simulations above. Again, we assume a telescope diam¬ 
eter of 39m, with an effective collecting area of 976 m 2 , and a 
total telescope+instrument throughput of 15%. Since we do not 
have a detailed simulation of the telescope PSF at optical wave¬ 
length, we approximate the theoretical, monochromatic PSF of 
the E-ELT at optical wavelengths by simply scaling the angular 
size of the diffraction-limited PSF at L-band. We do not consider 
a polychromatic PSF because this would increase the computa¬ 
tional time by an order of magnitude, and would only result in 
a low-frequency modulation of the observed spectra per pixel, 
which can easily be filtered out. The effect of seeing has been 
taken into account by scaling the amplitude of the PSF to 30% 
according to the Strehl ratio, and adding a seeing limited halo 
with a relative strength of 70%. This halo was modeled as a 
Moffatt proh 1(0 with /j = 2.5 and a = 265.4, to match a seeing 
FHWM of 0.6 arcseconds. 

As above, we constructed a three-dimensional data cube, 
with x-, and y-position, and wavelength on its axes. The spec¬ 
trum of the star is calculated, assuming 10 hours of observing 
time with the E-ELT, in photons per wavelength step, and trans¬ 
lated to photons per wavelength step per pixel using the normal¬ 
ized PSF. The same procedure is followed for the planet signal, 
except that it is first scaled according to the planet-star contrast 

3 Moffat profile is F(r) = F 0 [l + (r/a) 2 ]^ 
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(in this case 6.0xl0 -7 ), shifted in velocity by +30 km sec -1 , and 
added to the data cube at the calculated planet position. Since 
we assume that the starlight can be reduced down to the photon 
noise, to emulate noise only the square-root of the number of 
photons per pixel for the star is added to each pixel, and mul¬ 
tiplied by a random selection from a Gaussian distribution with 
cr = 1. As in the case of infrared HDS+HCI, we assume an ide¬ 
alistic case. Despite much success with HDS alone, the ability 
to remove stellar contributions to < 10 4 when combining with 
HCI remains to be proven because of the unknown behavior of 
stellar speckles at high-spectral resolution. Although some re¬ 
gions in the optical wavelength regime are affected by telluric 
absorption, for simplicity we ignore this for these optical simu¬ 
lations. Including it should not make a significant difference to 
our S/N calculations. The sky background is also ignored, since 
it is more than 23 magnitudes fainter (< 5 x 10 10 ) per pixel than 
the star at a good observing site. 

Each pixel location is subsequently cross-correlated with the 
stellar spectrum to produce a cross-correlation data cube, with 
x,y position and velocity (-150 to +150 km sec -1 ) on its axes. 
The slice at the planet velocity is shown in Figure 7. Since 
the cross-correlation technique automatically renormalizes the 
noise, the star is not visible in this data cube. The planet, at a 
contrast level of 6.0x10 7 at 25 milli-arcseconds from the star, 
is detected at S/N~10. 

6. Summary and outlook 

We have conducted simulations of HDS+HCI observations with 
the E-ELT, in the infrared and optical wavelength regime. The 
infrared observations target absorption lines in the thermal emis¬ 
sion spectrum of the planet. The METIS instrument will have the 
high-dispersion (R= 100,000) IFS capabilities to perform such 
observations, and we show that at M-band, rocky planets in the 
habitable zone of a Centauri A are within the realm of this in¬ 
strument. Since the angular separation of the habitable zone for 
this star is quite far out, it is expected that such planet will be 
in the sky background limited regime for the E-ELT. It may well 
be that more classical high-contrast imaging techniques, mak¬ 
ing use of angular and/or spectral differential imaging (ADI & 
SDI) will be similar, or even more powerful than HDS+HCI un¬ 
der these conditions since they make a more efficient use of the 
available photons. For more distant stars the relative contribu¬ 
tion of the sky background is even larger, but the planets will be 
closer in in angular space, towards the regime where ADI and 
SDI are less effective. In addition, HDS+HCI observation will 
provide the radial component of the planet orbital velocity, line 
strength of molecular absorption in the planet atmosphere, and 
can be used to determine the planet spin (see below). 

Optical HDS+HCI is technically more challenging, because 
of the more stringent requirements on the telescope adaptive op¬ 
tics, but it will be more powerful than infrared observations be¬ 
cause the sky background plays a significantly smaller role. We 
show that with a Strehl ratio of 0.3, one night of observation of 
Proxima Centauri with the E-ELT could reveal a rocky planet 
in its habitable zone, with a few dozen other late M-dwarf that 
could be surveyed in this way. 

The discovery and characterization of the planetary systems 
around our stellar neighbors constitute fundamental measure¬ 
ments of great scientific and philosophical value. The HDS+HCI 
technique has the potential to detect and characterize temper¬ 
ate rocky planets in their habitable zones. Exoplanet scientists 
should therefore not shy away from claiming a large fraction of 
the time on the future ELTs to make such observations. 
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First HDS+HCI observations 


HDS+HCI observations of unknown planetary systems will re¬ 
quire a high-dispersion (R~ 100,000) integral field spectrograph 
with adaptive optics capabilities. To our knowledge such a sys¬ 
tem does not exist on current telescopes. However, if the po¬ 
sition of the planet is known, then a slit spectrograph can also 
be used, orientated in such way that it encompasses both the 
planet and the host star. We tested the HDS+HCI technique 
on the well-known young planetary system /3 Pictoris b, using 
the R= 100,000 CRIRES spectrograph on the VLT, which has 
been presented in a separate publication (Snellen et al. 2014| >. 
Although the observations were conducted under mediocre (1- 
1.3 arcsec) seeing conditions, the starlight at the planet posi¬ 
tion was reduced by a factor of 8 to 30 compared to the peak 
of the stellar profile using the MACAO adaptive optics system. 
The 2.3 n m bandhead of carbon monoxide was targeted, and the 
planet, ~ 4000 times fainter than the star at this wavelength, was 
detected with an S/N of 6.4 in only one hour of observations 
including overheads. For comparison, r Bootis b, at a similar 
planet/star contrast and with the star of the same brightness as 
P Pictoris, was detected at a similar significance with classical 
HDS ( |Brogi et al.|2012[ > but in 18 hours of observing time. The 
factor of ~20 gain in observing time is as expected from the 
factor of 8-30 reduction in starlight at the planet position using 
HDS+HCI. 

As a testemony to the potential power of the HDS+HCI tech¬ 
nique, the observations revealed the planet orbital velocity, and 
its spin - a first for an extrasolar planet. 
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